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We consider a real massless scalar field in 3 + 1 dimensions satisfying a Robin 
boundary condition at a nonrelativistic moving mirror. Considering vacuum as the 
initial field state, we compute explicitly the number of particles created per unit 
frequency and per unit solid angle, exhibiting in this way the angular dependence of 
the spectral distribution. The well known cases of Dirichlet and Neumann boundary 
conditions may be rcobtaincd as particular cases from our results. Wc show that the 
particle creation rate can be considerably reduced (with respect to the Dirichlet and 
Neumann cases) for particular values of the Robin parameter. Our results extend 
for 3 + 1 dimensions previous results found in the literature for 1 + 1 dimensions. 
Further, we also show that this inhibition of the dynamical Casimir effect occurs for 
different angles of particle emission. 

PACS numbers: 03.70.+k, ll.lO.-z 
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I. INTRODUCTION 

The dynamical Casimir effect (DCE) basically consists of the emission of quanta from 
a moving body in vacuum due to its interaction with a quantized field |i"5j. Another 
manifestation of the DCE, which is a direct consequence of the particle creation phenomenon 
if we invoke the energy conservation law, is a radiation reaction force acting on the moving 
body. This dissipative force gives rise to an irreversible exchange of energy between the 
moving body and the quantized field. In other words, the energy dissipated from the moving 
body is converted into real excitations of the quantized field, i. e., real particles. One can also 
understand the DCE in the opposite way, namely, from the fluctuation-dissipation theorem 
[HI Ej and from the fact that the static Casimir force acting on a fixed body, though zero, 
has non- vanishing fluctuations [H]. In this scenario, one expects that a moving plate may 
be acted by a dissipative force (under certain circumstances) which is proportional to the 
fluctuations of the Casimir force on the static plate [QHTT] (for the case of a moving sphere 
see Ref. [12]). 

However, the quantized field and the moving body can also exchange energy reversibly, 
which means that the force exerted on the moving body acquires in this case a dispersive 
part, as it occurs when Robin boundary conditions (BC) are considered |13l HI]. During 
the first two decades after the pioneering paper by Moore [IJ, the calculations on the DCE 
were usually done with scalar fieds. The consideration of electromagnetic fields was made 
by the first time in 1994 [15] (see also Ref(s) [161 IE]); ^ind since then great attention has 
been devoted to the DCE. Detailed reviews on the DCE can be found in Refs. [HI [T9] . 
Several experimental proposals to observe the DCE have been made in the last years. We 
shall briefly comment on a couple of them (for more details see, for instance, Dodonov's 
paper [20]). 

The so called motion induced radiation (MIR) experiment [21] is based on the simulation 
of a mirror's motion by changing the reflectivity of a semi-conductor by irradiating it with 
appropriate laser pulses, an ingenious idea firstly introduced by Yablonovitch in 1989 [22] in 
a paper where the main concern was to propose ways of simulating highly accelerated frames 
in order to enhance the Unruh radiation. A few years later, this same idea of creating a 
dense electron-hole plasma in a thin semiconductor by irradiating it with laser pulses, was 
also discussed by Lozovik et al [23]. Though there are many promissing aspects in the 
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MIR experiment, it is worth mentioning that the MIR experimentalists may have to deal 
with some difficulties. A ffist one is related to the limitations in the signal-to-noise ratio 
present in their experiment caused by thermal effects, in case they run the experiment at 
4.6 K, as pointed out by Kim et al [21]. A second one is the inffuence of damping in a 
parametric amplification process. And in the MIR experiment, the electric permittivity of 
the semiconductor slab after excitation by the laser pulse acquires a non-negligible imaginary 
part [25j, so that dissipation effects are inevitable. In this shown by Dodonov, 

damping plays an important role and the emergence of a superchaotic quantum state may 
occur leading to a highly superPoissonian statistics for the distribution function of quanta 
[26]. 

Another interesting proposal was made by Kim et al [27J. They suggest that an indi- 
rect measurement of the dynamical Casimir photons generated by means of the mechanical 
motion of a film bulk acoustic resonator is detected with the aid of a superadiance mecha- 
nism. More recently, Dezael and Lambrecht [28] proposed that a (dynamical) Casimir-like 
radiation may arise from an effective motion of mirrors obtained by the interactions of an 
optical parametric oscillator with a thin non-linear crystal slab inside. In 2011, Kawakubo 
and Yamamoto [29j proposed that photons could be generated by means of a non-stationary 
plasma mirror, being the photons detectable by an excitation process of Rydberg atoms 
through the atom-field interaction. Still in 2011, Faccio and Carusotto [SD] proposed a pho- 
ton generation mechanism in the near-infrared domain obtained by a train of laser pulses 
applied perpendicularly to a cavity, made of non-linear optical fiber, which modulates in 
time the refractive index of the medium filling the cavity. 

Finally, forty years after its theoretical prediction made by Moore [1], the ffist experi- 
mental observation of the DCE was announced by Wilson and collaborators [31], in an ex- 
periment where these authors use a superconducting circuit consisting of a unidimensional 
coplanar transmission line with a tunable electrical length. The change of the electrical 
length is performed by modulating the inductance of a superconducting quantum interfer- 
ence device fixed at one end of the transmission line. This modulation is achieved with the 
aid of a time-dependent magnetic flux through the superconducting quantum interference 
device. The electromagnetic field along the transmission line is described in terms of a field 
operator given by a scalar field 4>{t,x) obeying a massless Klein-Gordon equation in 1-1-1 
dimensions and submitted to a Robin BC with a time-dependent Robin parameter ^[t) in 
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the following manner [32] 



0(t,O)^7W(5.</')(t,O). 



(1) 



This kind of BC was also considered in Ref(s) [331IM]. In the context of the DCE, Robin BC 
appeared for the first time in the papers by Mintz and collaborators [I3l[35], who investigated 
a real massless scalar field (j){t, x) satisfying a Robin BC at a moving plate when observed 
from an inertial frame in which the plate is instantaneously at rest (we shall refer to this 
frame as tangential frame), according to the formula 



where the prime superscript is to remind us that the BC is taken in a tangential frame and 
7 is now a time-independent Robin parameter. 

For an oscillating mirror that imposes Dirichlet BC on a scalar field in 1 + 1 dimensions, 
the total particle creation rate is a monotonic function of the mechanical frequency of the 
mirror [36]. However, it was shown in [35] that this is not always the case when the field 
obeys a Robin BC given by Eq. Q. For a given value of 7, there is an interval in which an 
increase in the oscillation frequency of the mirror leads to a decrease in the particle emission. 
This can be understood as a kind of "decoupling" between the mirror and some of the field 
modes. 

One important question that remains is whether or not this interesting effect still occurs 
in 3 + 1 dimensions. This is not a trivial issue, once the phase space available for the field 
in (i = 3 + 1 is much larger than that in 1 + 1 dimensions. To answer this question is the 
main purpose of our paper and we shall do that by considering a real massless scalar field in 
3 + 1 dimensions satisfying a Robin BC at a non-relativistic moving mirror (in the tangential 
frame). As far as we know, all papers about DCE with Robin BC deal with models in 1+1 
dimensions, except one [37], which did not discuss the problem we are interested here. 

An extra motivation to study this model is the connection between the DCE for a real 
massless scalar field in 3 + 1 dimensions and the DCE for the electromagnetic field inter- 
acting with a perfectly conducting plate. The latter problem can be separated into two 
problems: a vector potential representing the transverse electric (TE) polarization, which 
is associated to a Dirichlet BC, and a vector potential representing the transverse magnetic 
(TM) polarization, which is associated to a Neumann BC. Since the parameter 7 allows a 
continuous interpolation between Dirichlet and Neumann BC, we expect that in the limit 



(2) 
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7 — !■ the DCE for the massless scalar field coincides with the TE polarization contribu- 
tionto the electromagnetic DCE, while for 7 — )■ oo the TM polarization contribution to the 
electromagnetic DCE is recovered. 

The structure of this paper is as follows. In Sec. II, we calculate the relations between 
the creation and annihilation field operators in the remote past ( "in" operators) and in the 
far future ("out" operators), in the Heisenberg picture. In Sec. Ill, we show our results for 
the particle emission rate for a specific but typical motion of the mirror. Finally, in Sec. IV, 
we discuss our results and make a few comments on their possible consequences. 



II. BOGOLIUBOV TRANSFORMATION AND PARTICLE SPECTRUM 

Let us consider - in the Heisenberg picture - a massless scalar field in 3 + 1 dimensions 
written in terms of the time-dependent operators a^(t,k) and a(t,k), where k is the wave 
vector. In the distant past, these operators are relabeled as the creation and annihilation 
operators cij„(k) and aj„(k), whereas in the far future they are relabeled as alu^{k) and 
aouf(k). The time evolution of the operators a'' and a depends on the interaction between 
the field with an external agent (modeled, in the present paper, by a moving boundary). 
The "out" operators can be expressed as a combination of the "in" operators via Bogoliubov 
transformations [38] 

ao„i(k) = a(k)ai„(k) + /3(k)aJ„(k), (3) 

where a and /3 are named Bogoliubov coefficients. Assuming that in the remote past the 
system is in the vacuum state |0), the spectral density of the created particles after the 
movement of the mirror has ceased is given by 

|^(k) = ^(0|aL(k)a.„,(k)|0) = |/3(k)|l (4) 

Notice that there will be particle production if, and only if, /3(k) 7^ 0, i.e., if the annihilation 
operator aoui(k) is "contaminated" by the creation operator aj„(k). The relation between 
the "in" and "out" operators can be calculated exactly for the case 1 + 1 dimensions [5J with 
Dirichlet or Neumann EC. However, only approximate approaches are currently known for 
higher dimensional space-times. If the movement of the mirror is non-relativistic and has a 
small amplitude, the perturbative method introduced by Ford and Vilenkin pQ] applies. In 
the following, we will use the Ford- Vilenkin approach to find the Bogoliubov transformation 



6 



for the massless scalar field satisfying the wave equation 



(5) 



and obeying a Robin BC imposed by a moving mirror. We start writing the field operator 
as the field under static BC plus a perturbation, that is 

0(t,r) = 0o(t,r) + 50(t,r), (6) 

where we assume that the undisturbed field 0o obeys the static boundary condition BC 

00 (t,r)- 7^.00 (t,r)l =0, (7) 

- 2 = 

where 7 is a time-independent parameter. 

The perturbation 50 (which is assumed to be small) gives the first order contribution to 
the total field 0(t, r) caused by the movement of the mirror and therefore will be responsible 
for the emergence of the DCE. 

Let z = 5q{t) be the position of the mirror at a given instant. We assume that, in the lab 
reference frame, the mirror starts at rest at z = 0, undergoes a given prescribed movement 
and finally settles down at z = for large times. The perturbation 50 will be small as long as 
the speed of the mirror is continuous with continuous derivatives and is much smaller than 
the speed of light, |5g'(t)| ^ 1, in natural units. Once we consider a bounded movement for 
the mirror, we also require that it possesses a small amplitude. More specifically, this means 
that \5q\ ^ 1/wo) where Uq is the dominant mechanical frequency and this assumption will 
allow us to neglect terms 0[{5q)% C[(50)^], O[505g] and 0[{5q)% 

We shall follow the same procedure as in [35j. At the tangential frame at a given instant, 
the BC is given by 

0'(t',r')=7'|^(t',r'), (8) 

where 7' is a t'-independent parameter. This BC can be cast in the lab frame with the help 
of the appropriate Lorentz transformations. For non-relativistic velocities, one may expand 
in first order in 5q{t) to find 

9,0 (t, r) + 5g(t)9t0 (t, r) - 7- V (t, r)l =0, (9) 

- z=5q{t) 

where, in this approximation, we consider that the Robin parameter is not affected by the 
Lorentz transformation, so that 7' = 7. 
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We substitute (|6]) in ([9]) and expand the result thus obtained in the small parameters 
retaining only terms up to linear order in 5q and 5(f), finding 

-6q{t)dt<Poit,r)U=o- (10) 

Therefore, the perturbation 6(j){t, r) obeys a time-dependent BC at a fixed position {z = 0), 
which is associated to the motion of the mirror via Sq{t). 

It is now convenient that we express the field in the frequency domain using a Fourier 
transform, as follows. Let us define $(a;, ky; z) as the Fourier transform of the field r), 
so that 



'2„ Aujt „—iU.\i-ri 



dt I d\e"^'e 



0(t,r). 



(11) 



Notice that $(a;,k||; z) obeys the Helmholtz equation 



{u^-kl + d',)<l>{u,kf,z)=0. (12) 

Once the massless free scalar field is a solution of the d'Alembert equation (|5|, 

= [{u + tef - A;J] , with e ^ 0+. (13) 

With this definition, is a complex function of u, with a branch cut along the real axis 
between — fcy < a; < fcy [16j, where k^\ = |k|||. 

The undisturbed field is the solution of the wave equation ([s]) subject to the static BC 

0: 



0o(t,r) 



/o V 4:71-^ |wk| yl + 7^5 



sm{kzz) + '-ykz cos{kzz) 



(14) 



where the field normalization is chosen so that the creation and annihilation operators obey 
the commutation relations 



[a(k),at(k)'] = (27r)3(5(k-k'). 



(15) 



We can now use (11) to write down the Fourier transform of the undisturbed field. 



$o(w,k||;z) = y4(u;, fc2;7) sinik^z) + ^k^cosik^z) a(k)9(a;) - a^(-k)6(-a;) 6(fc^),(16) 
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with the normahzation factor 



(17) 



The Fourier transformations of the perturbation r) and the mirror's law of motion 
5q{t) are written as 



5$(a;, k|j; 2;) = / dt j d 
J —00 J 



'2^^^gii^tg-iki|.r 



ii-"ii50(t,r) 



and 



dte"^'6q{t). 



Notice that 5$(tj,k||; z) obeys the Helmhohz equation 



(18) 



(19) 



(20) 



Recalhng that we must take only the solution that propagates outwards from the moving 
mirror, the solution of the previous equation for z > is given by 

1 



ik 



-dzS^{cu, k||; 0)e 



ikzZ 



(21) 



We now notice also that 6^{uj, ky; z) obeys the BC given by the Fourier transform of (10), 
namely, 

dco' 



2n 



a,(5$(a;,kjj;0) -7-^5$(w,k|j;0) = 7"^ 

^■yuu' - 7kjj j $0 (w k|| ; 0) 
Once the Fourier transform of the Ford-Vilenkin ansatz ^ is 

$(w,k||;2;) = <l>Q{u,ki\;z) + 6<^{uj,kf,z), 



SQ{cu-uj') d,^o{uj',kfO)+ 



(22) 



(23) 



we now need to find an expression for $(u;, ky; z) that allows us to relate the Fourier trans- 
forms of the field in the remote past and distant future. This can be done straightforwardly 
by using Green's functions, as in [35]. 

Let us start with the one-dimensional version of Green's identity, dz{g[dzf] — f[dzg]} = 
g[d'^f] — f[d],g]. We now identify the function / in this formula as the perturbation 5$ and 
the function g as the Green's function of the Helmholtz operator, i.e., a function such that 
(^tu^ — fcy + d^^ g{uj, k||; z, z') = S{z — z'). After integration by parts, we can show that 



5$(u;, k|| ; z') = g{^. k|| ; 0, z')dM{.uj. k,, ; Q)dzg{uj, k,, ; 0, z')5^uj, k,, ; 0) 



(24) 
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We may freely define the Green's function g so that it obeys the Robin BC at the surface 
of the static mirror, i.e., g{uj, ky; 0, z') = 79^(7(0;, ky; 0, z'). Then, the identity (24) leads to 

^(w, k|| ; z')g{uj, k,, ; 0, z')\dM{uj, k,| ; 0) - 7"^5$(u;, k,, ; 0)1 (25) 



and, using Eq. (22), we can write the perturbation as 

""^ du' 



6^UJ,^;z')=^-'g{uj,kf,0,z') 



2n 



6Q{uj-uj') 9,$o(c^',k||;0) + 



'yuu' — 7kn )$o(c(;', ky; 0) 



(26) 



Replacing (26) in (23) one finds 



1 



$(a;, k||; z') = $0(0;, ky; 2;') + -g{u, ky; 0, z' 

II n ^ H 

(^-fuju - 7kyj<l>o(a;',k||;0) 



dcu' 
27 



(27) 



The field $(u;, ky; 2;') can be written in two different forms. The first one involves $m, 
the Fourier transform of the unperturbed field (/)o(t, ry; 2;') in the distant past, that is, the 
"in" field, to which we associate retarded Green's functions. Then, 



<l>(cu,ky;z') = $in(c<;,ky; z') + -g,^ti(^,kn;0, z' 



7 



27 



7ww' - 7kM )$o(i^', ky; 0) 



(28) 



where $in(c<;, ky; 2;') is the Fourier transform of 0in(t,r), which is simply the unperturbed 
field at the remote past, 

(29) 



0o(t,r) ~ 0m(t,r). 

t— >— 00 



Analogously, one can use advanced Green's functions to express the field $0(1^! ky; z') using 
the "out" field configuration, as follows. 



^{uj,kf,z') = <l>out(^,ky; z') + -g^dv{oj,kfO,z') 
7ww'-7kj)$o(a;',ky;0) 
where ^out{^, ky; z') is the Fourier transform of 



4>o(t,r) ~ 0out(t,r). 



27 



5Qiu-uj') 9,$o(a;',ky;0) + 



(30) 



(31) 



From Eqs. (30) and (28), we can relate the "in" and "out" field operators through the 
advanced and retarded Green's functions. 



$out(^ky; 2;') = $in(w,ky;2') +- ^iret (w, ky ; 0, 2;') - fi'adv(w, ky ; 0, z') 



X 
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— 5Q(u; - uj')\d,^o{uj\ k||; 0) + (^iujuj' - 7kJ j %{uj\ kj,; 0)J (32) 



Using the solution (16) of the Helmhohz equation for the field subject to Robin BC, one 



can reexpress (32) as 

$out(c^,k||;2;') = $in(w,k|j;2;') + 



2i 



1 + -i^kl 

^Aico', ^)6Q{co - co')k'z (l - 7^k5 + 7 W) x 
a(q)e(fc;)-at(-q)e(-A;;)le(fc: 



where we defined 



with 



q := k|| + k'J, 



211/2 



(33) 



(34) 



(35) 



Once again, using (16) it is straighforward to show that 

2ik, 



aout(k) = ain(k) + 



~ duo" 



\U0" 



v/|c^|(l + 7'A;2) 27r \/ (1 + 7^^^) 

aM'mK) - «L(-q")0(-A::)l0(A:f ), 



6Q{uj - uj") (l - 72kJ + -i^ujJ'^ X 

(36) 



with q" = kj| + k'!,z. The last expression is a linear relation between the annihilation field 
operator in the far future with the creation and annihilation operators at the remote past. 



That is, (36) is the Bogoliubov transformation we were looking for. 



As advertised, the vacuum state is annihilated by the "in" operator ain(k), but not by 
the "out" operator ao„t(k) due to the presence of a]^„(k) in the r.h.s. of (36). This clearly 



indicates a nonzero particle creation in 3 + 1 dimensions due to the DCE for a massless field 
subject to Robin BC at a moving mirror. 

After some straightforward manipulations, from Eq. (|4]) and using the Bogoliubov trans- 
formation ([36|), we can finally calculate the spectral density of created particles: 



dN{\i) 



A 



^kl 



dcu' \6Q{co + 
271 ' l+lokf 



(l-fkJ-7W)', (37) 



d^k (27r)^|u;|(l+7ofc2)^|k„| 
where A is the area of the moving plate. Taking into account the axial symmetry of the 



problem, it is interesting to rewrite (37) in terms of the variables 



(38) 
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with = |k| and 9 G [0,7r/2]. Analogously, we have 

k', = Jlo'^ - kl = Vuj'^ - oj^ sin^ 9, 



(39) 



which leads us to our expression for the particle spectrum per unit area, per unit solid angle 



1 dN 



1 



Au'^ cos^ 9 



A dudVL {2'Ky 1 + i^uj"^ cos2 9 J^.^^g 2ti 



(^1 - 7 V sin^ ^ - 7W) . 



(40) 



The expression (40) is valid for an arbitrary motion of the mirror, as long as it is sufficiently 



slow {\6q\ <^ 1) and has a small amplitude. In the following section, we present our results 
for a specific but very useful law of motion, that of a mirror which oscillates, practically, at 
a single frequency. 



III. RESULTS AND DISCUSSION 



Let us consider that the movement of the mirror is described by 



5q{t) = eocos(wot)e"'*'/^. 



(41) 



which is a typical motion considered in investigations of the DCE [16], where eo represents 
its amplitude, Uq the dominating mechanical frequency and r the effective time interval of 



the oscillation. The Fourier transform of Eq. (41) is 

1 



SQ{u) = eoT 



+ 



1 



In the limit with uqt ^ 1, we obtain 



\6Q{u 



|2 _ ^.2 



ot[5{uj - uo) + 6{u + cjo)], 



(42) 



(43) 



which is reasonable once (42) possesses two narrow peaks around tu = ±ujo in this limit. 



Substituting (43) into (40) and defining Af the normalized particle spectrum per unit area, 



per unit time, per unit solid angle, given by 
we can write 



1 dN felr 



AdudVt VStt 



M{(jj, cjo, 7, = ^0 -^("i 



(44) 



(45) 
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where 

0[1 -a(l + sin^)], (46) 



with a = io/uJo and /3 = uJo^y. Notice that the Heaviside function present in Eq. (45) 
indicates the following interesting features of the particle spectrum: particles associated to 
a certain frequency u can be created in all angles if < a; < a;o/2; for particles associated 
to a frequency u > ijJo/2, there is no particle emission for angles larger than 9o{u), with 

9o{uj) = arcsin (^^^^^) • (47) 



Moreover, observe that expression (45) coincides with the results found in [TB] in two partic- 



ular cases. In the limit 7 — )■ 0, Eq. (45) reproduces the particle spectrum for TE photons, 
while for 7 — > 00 it corresponds to the spectrum of TM photons. 

The normalized particle spectrum per unit area, per unit time, here labeled as A/", can 



be obtained by integrating (45) in the solid angle, with angular range 9 E [0, 7r/2]: 

Ar(u;,u;o,7)=^o-^(«,/9) (48) 

with 

J^{a,f3) = 27i d9 sm9 :F{a,f3,9). (49) 
Jo 

The results for uq = 1 and different values of 7 in f/{u,uo,'~f) are shown in Fig. 1. The 
spectrum is symmetric around u)q/2 for every value of 7. This is in agreement with the 
expectation that the particles are created in pairs with frequencies Ui and U2 so that ui+ijJ2 = 
Uq, the driving frequency of the mirror fTOJ. 

One can also analyze the angular spectrum of emitted particles. Our results for N", which 



follow from Eq. (45) for a particular frequency u = uJo/2 and = 1, are shown in Fig. 2. 
A few comments are in order. Firstly, for small values of the Robin parameter, 7^1, the 
emission is mainly forward, with large angle emission being suppressed. In second place, for 
7 ~ 2, the emission is strongly inhibited for every angle, but it also acquires a maximum 
around 9 ~ Irad. If the 7 is further increased, the emission rate rises again as a whole, but 
the particle production is highest around 9 ^ Irad, and not around the normal direction 
9 = 0. Another interesting and subtle point becomes evident when we compare the angular 
spectrum in the case of Robin boundary conditions with 7^1 and that of the Neumann 
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FIG. 1. (Color online.) The normalized spectrum J\f in arbitrary units, as function of oj, for 
different values of 7, with ujq = 1. The particle emission at a given frequency u is initially reduced 
as 7 increases from zero until a given angle and then it starts increasing, reaching a fixed plateau 
for 7 — 00. 



BCs. We see from Fig. 2 that there is no emission for — )■ 7r/2 no matter how large the 
values of 7 are. This is basically due to the factor cos^ 6 in (45). However, as shown in 
[16] . there is a finite emission rate at ^ — t- 7r/2 for Neumann BC. The apparent paradox can 



be solved by noticing that the limits 6 ^ n /2 and 7 —i- 00 in Eq. (40) do not commute. 
Fortunately, this subtlety affects only grazing angles and the limit 7 — ?■ 00 can be identified 
with the Neumann BC whenever 6' = 7r/2 is not the only angle in consideration. 

Finally, let us turn our attention to the total number of particles emitted N, which is 
given by 

Niuo,-f)=ulF{f3) (50) 

where 

POO 

F{(3) = / daF{a,(3). (51) 

The ratio between the emission rate Nr = N{uo,'j) and the emission rate = N{uo,0) 
(Dirichlet BC), is given by 

NR/Nn = F{f3)/F{0), (52) 
which is showed in Fig. 3. As first demonstrated in [16], the total emission rate for Neumann 
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Polar emission angle [rad] 



FIG. 2. (Color online.) Angular spectrum N of emitted particles with frequency oj = ooq/2, with 
cjQ = 1, as a function of the emission angle (measured with respect to the normal to the mirror) 
for various values of 7. Solid (black) line: 7 = 0; dashed (red) line: 7 = 1; long dashed (green) line 
: 7 = 2; dotted (blue) line: 7 = 3; dashed-dotted (magenta) line: 7 = 6. 



BC is 11 times larger than that for Dirichlet BC in 3 + 1 dimensions. Our results not only 
confirm the factor 11 but also show that Dirichlet and Neumann cases are connected by a 
non-monotonic curve. Indeed, the Robin BC with 7CJ0 — 2 provides a very strong inhibition 
of the particle production. This property of the DCE with Robin BC had already been 
noticed in 1 + 1 dimensions [35], but it was not obvious a priori that this property would 
be preserved at higher dimensions. 



IV. FINAL REMARKS 

We have investigated the DCE for a real massless scalar field in 3 + 1 dimensions satisfying 
a Robin BC at a non-relativistic moving plate (in the tangential frame). We used the 
perturbative method of Ford and Vilenkin [32], valid for non-relativistic motions of small 
amplitudes to evaluated the Bogoliubov transformations between creation and annihilation 
operators in the remote past, aj„, and distant future, aout- This allowed us to determine 
the spectral and angular distributions of the created particles caused by the moving mirror. 



Eq. ( 40 ) exhibits our expression for the particle spectrum per unit area, valid for a general 



15 
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FIG. 3. (Color online.) Total emission rate for Robin BC normalized by the Dirichlet rate. Notice 
that the Neumann limit (7 — )• 00) is correctly recovered. The strong decoupling between the mirror 
and the quantum vacuum field at 7^0 ^ 2 is typical of the Robin BC [T3l [35] . 



(non-relativistic) law of motion for the moving plate. Assuming the oscillating law of motion 
given by Eq. (liTl) - a typical motion considered in investigations of the DCE - we obtained 



an explicit expression for the spectrum per unit area, given by Eq. (45). In the limit 7 — )■ 
(Dirichlet case) we recovered the particle spectrum for TE photons, whereas for 7 — )■ cxo 
(Neumann case) we reobtained the spectrum of emitted TM photons. 

Our results also show that, although in the limits 7 — > and 7 — )■ 00 the total number 
of created particles is a monotonic function of the mechanical frequency of the plate, Uq, 
the same is not true for intermediate values of 7. In fact, for any fixed positive 7, the 
total number of created particles is not a monotonic function of uo- More than that, the 
strong inhibition of the DCE that occurs in 1 + 1 dimensions, when the Robin BC is present, 
also occurs in 3 + 1 dimensions. Indeed, the total number of created particles shown in 
Fig. 3 is dramatically reduced for a quasi-harmonic motion with a frequency Uq such that 
71^0 — 2, where 7 is the parameter that characterizes the Robin condition. Naively thinking, 
this surprising effective decoupling between the plate and the quantized field, predicted 
previously in 1 -|- 1 dimensions [IHl [35] was not expected in 3 -|- 1, since in the latter case 
only the field modes that propagates perpendicularly to the plate are expected to behave as 
the field modes in 1 -|- 1 dimensions. 
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The suppression just discussed show how the dynamical Casimir effect may be strongly 
dependent on the boundary conditions employed. This fact is very important, since any 
extra information about the created particles, either in the total number of them or in 
the angular dependence of the corresponding spectral distribution, can be extremely useful 
to better identify the dynamical Casimir photons in a given experiment. Since the Robin 
parameter 7 can be interpreted as the plasma wavelength of a given material, the results 
presented here may be of some help for future experiments, at least as a source of concern 
and caution when dimensioning these experiments. For the moment, realistic values for the 
product juq are still very far from a strong suppression, but this may not be the case in 
future experiments. It would be interesting if the peculiar signature of the DCE with Robin 
conditions - the strong inhibition in the particle emission - could eventually be captured 
in experiments. Considering the connections between Robin BC and the theoretical model 
underlying the first experimental observation of the DCE [31], we believe that a thorough 
study of the implications of the Robin boundary conditions on the DCE is crucial. 

As a final comment, since the dissipative force acting on a moving mirror is closely 
related to the particle creation, we expect that this dissipative force will also suffer a similar 
inhibition in 3 + 1 dimensions, as it occurs in the 1 + 1 dimensional case [13]. This calculation 
will be left for a future work. 
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